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的 调节 (mediated moderation, meMO) 模 型 关于 误差 方差 齐 性 的 假设 经 常 被 违背 ， 
中 也 缺乏 测量 meMO 效应 大 小 的 指标 。 对 于 单 层 数据 ， 本 文 借助 于 两 层 建 模 的 思想 ， 
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应 用 研究 
提出 了 一 种 可 用 于 处 理 方差 


非 齐 性 的 两 层 有 中 介 的 调节 (2meMO) 模 型 ; 给 出 了 
节 效 应 大 小 的 效应 量 。 


于 测量 meMO 分 析 中 总 调节 效应 、 直 接 调 节 效 应 和 有 中 介 调 
通过 Monte Carlo 模拟 研究 ， 比 较 了 meMO 和 2meMO 模型 在 参数 和 效应 量 估计 上 的 表现 。 


并 通过 实际 案例 解释 了 2meMO 模型 的 应 用 以 及 效应 量 的 计算 和 解释 。 
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引言 


在 很 多 领域 的 实际 研究 中 ， 中 介 (mediation) 效 
应 分 析 和 调节 (moderation) 效 应 分 析 是 探讨 自 变 量 
了 和 因 变 量 了 之 间 复 杂 关 系 常用 的 两 类 方法 。 随 着 
人 研究 问题 的 深入 ,在 更 系统 的 理论 框架 下 , 将 中 介 
效应 与 调节 效应 结合 ,探讨 变量 之 间 更 为 复杂 关系 
的 研究 方法 越 来 越 受到 重视 ( Kwan & Chan, 2018; 
叶 宝 娟 , 温 忠 诅 ，2013)。 其 中 ， 有 中 介 的 调节 
(mediated moderation, 简称 meMO) 模 型 在 探讨 调 
节 变 量 作 用 机 制 方面 , 提供 了 强 有 力 的 分 析 方 法 。 

meMO 模型 分 析 的 重点 是 探讨 调节 变量 Z 对 自 
变量 (0) 与 因 变 量 () 之 间 关 系 的 作用 机 制 ， 如果 Z 
对 和 和 了 之 间 关 系 的 调节 通过 中 介 的 调节 变量 
(M)(mediating moderator) 起 作用 ， 则 称 Z 通 过 MM 对 
X Fil Y 之 间 关 系 间接 的 调节 为 有 中 介 的 调节 
(meMO)。 基于 不 同 理论 的 假设 , meMO 模型 近年 来 
在 心理 学 研究 中 得 到 了 广泛 的 应 用 (如 ，Baron & 
Kenny, 1986; Liu et al., 2012; Muller et al., 2005; Œ 
玲 晓 ^E. 2019; 杨 文 圣 SF, 2019; Hie ^m, 


2020). 在 方法 研究 领域 ， 也 有 许多 研究 者 就 meMO 
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模型 的 分 析 步 又 和 检验 方法 进行 了 较为 详尽 的 介 
绍 和 总 结 (Hayes, 2018; Kwan & Chan, 2018; Ng et 
al., 2019; 叶 宝 娟 , ii EB, 2013)。 然 而 , 在 实际 应 
用 中 meMO 模型 仍 至 少 面临 以 下 三 个 方面 的 问题 。 

首先 ， 对 于 meMO 模型 ， 由 于 缺乏 合理 的 
meMO 效应 量 的 指标 ,应 用 研究 中 很 少 报告 其 效应 
量 ， 或 者 研究 中 所 选用 的 替代 性 指标 与 meMO BT 
究 感 兴趣 的 问题 不 匹配 。 针 对 研究 报告 中 过 分 强调 
统计 显著 性 的 问题 ， 越 来 越 多 的 期 刊 要 求 研 究 者 提 
供 效应 量 和 区 间 佑 计 的 结果 (如 , American Education 
Research Association, 2006; Cumming, 2014; Funder 
et al., 2013; Rozeboom, 1960; Wilkinson & American 
Psychological Association Task Force on Statistical 
Inference, 1999), 国内 的 《心理 学 报 》 和 《心理 科 
学 ) 等 杂志 也 有 同样 的 要 求 。 然 而 , 到 目前 为 止 , DF 
究 者 没有 结合 meMO 问题 的 特点 对 其 效应 量 进行 
针对 性 的 测量 和 解释 。 有 研究 者 借鉴 其 他 类 似 的 效 
应 量 作为 蔡 代 品 ， 如 采用 回归 分 析 中 常用 的 AR? fi 
述 增加 了 有 中 介 的 调节 路 径 后 ， 因 变量 变异 解释 率 
的 增加 (Aiken et al., 1991), 但 这 些 指标 并 不 合适 。 
因为 从 效应 量 定义 的 本 质 来 看 , 效应 量 是 反应 某 种 
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现象 (如 影响 程度 ) 大 小 的 数量 化 表示 , 合适 的 效应 
量 应 直接 针对 所 感 兴趣 的 研究 问题 (Kelley & 
Preacher, 2012)。 因 此 ,或 待 解决 的 问题 是 如 何 围绕 
meMO 分 析 感 兴趣 的 “Z 通过 M 调节 并 和 了 之 间 的 
关系 ”问题 ， 构 建 合理 的 效应 量 指标 ， 对 meMO 分 
析 中 关注 的 有 中 介 的 调节 效应 大 小 进行 评价 。 

其 次 ,有 中 介 的 调节 (meMO) 和 有 调节 的 中 介 
(moderated mediation, 简称 moME) 虽 然 解释 的 侧 
重点 不 同 , 但 是 在 实际 应 用 中 经 常 被 研究 者 混淆 。 
Baron 和 Kenny (1986) 最 早 提出 了 meMO 和 moME 
的 概念 ,很 多 研究 随后 就 这 两 类 分 析 进 行 了 详细 的 
前 述 (Edwards & Lambert, 2007; Liu et al., 2012; 
Muller et al., 2005). Hayes (2018) 结 合 概念 模型 和 统 
计 模 型 ， 对 meMO 和 moME 进行 了 详细 的 分 类 ,并 
提供 了 PROCESS 分 析 程 序 。 然 而 ,对 于 应 用 研究 
者 来 讲 清晰 区 分 两 类 模型 并 不 容易 , 人 研究 中 经 常会 
出 现 理论 假设 中 变量 之 间 的 关系 是 moME( 或 
meMO) 的 叙述 ， 而 在 统计 模型 结果 的 解释 上 却 采 
用 了 meMO( 或 moME) 的 解释 (如 ，Gonzalez-Mulé 
et al., 2016; Liang & Chi, 2013; Park & Searcy, 2012; 
Tang, 2016; Tang & Naumann, 2016; Zacher et al., 
2012; 杨 文 圣 等 , 2019)。 究 其 根本 原因 是 因为 在 数 
据 分 析 上 , meMO 与 moME 对 应 的 统计 模型 相同 ， 
导致 研究 者 对 meMO 和 moME 产生 了 混淆 。 具 体 
KDE, 常用 的 meMO 模型 分 为 meMO-I 型 和 
meMO-II 型 两 类 , meMO-I 模型 与 第 一 阶段 moME 
模型 等 价 , meMO-II 模型 与 第 二 阶段 moME 模型 等 
价 (Kwan & Chan, 2018)。 本 文 关注 的 问题 是 如 何 从 
概念 模型 出 发 ,在 统计 模型 的 表述 上 突出 memo 
模型 关注 的 变量 间 关 系 以 及 调节 变量 的 作用 机 制 ， 
定义 合理 的 效应 量 ,， 以 在 结果 解释 上 区 分 meMO 
和 moME 两 类 模型 。 

最 后 , meMO 模型 假设 误差 方差 齐 性 ， 回 归 系 
数 显著 性 的 检验 也 基于 这 一 假设 ， 然而 实际 应 用 中 
这 一 假设 常常 不 成 立 ,， 尤其 是 包含 调节 效应 的 模型 
(Aguinis & Pierce, 1998; DeShon & Alexander, 
1996)。 例 如 ,对 于 含有 交互 作用 的 回归 模型 

Y, = bytb X; bZ; +X; Z te, (1) 
式 () 中 Y, 和 和 2 分 别 表 示 因 变量 、 自 变量 和 调节 
变量 ; 表示 非 标准 化 的 回归 系数 ; YE 28 8; - 
N(0,0?), i =1, 2,…, n, n 表示 样本 量 。 方 差 齐 性 的 
假设 条 件 , 意味 着 对 于 自 变量 取 不 同 的 值 x; 时 ， 
误差 方差 相等 。 然 而 , 含有 显著 交互 作用 的 模型 很 
难 满足 方差 齐 性 的 前 提 假 设 ， 由 此 会 导致 第 二 类 错 


误 的 增 大 和 检验 力 的 下 降 ， 当 交互 作用 项 对 应 的 回 
归 系 数 显 著 不 为 零 时 ,其 参数 估计 对 应 的 标准 误 有 
偏 ， 从 而 导致 推断 统计 的 结果 (包括 效应 量 ) 出 现 偏 
差 (Aguinis & Pierce, 1998; DeShon & Alexander, 
1996)。 本 研究 拟 对 传统 的 meMO 模型 进行 拓展 ,在 
此 基础 上 定义 更 合理 的 测量 meMO 效应 量 的 指标 ， 
同时 也 解决 误差 依赖 于 自 变 量 的 方差 非 齐 性 问题 。 

为 了 处 理 调节 效应 分 析 容 易 导 致 的 方差 非 齐 
性 的 问题 ,Yuan 等 人 (2014) 采 用 两 层 模型 扩展 了 传 
统 的 调节 效应 回归 模型 ， 并 据 此 定义 了 测量 调节 效 
应 大 小 的 指标 。 虽 然 Yuan 等 人 (2014) 没 有 涉及 
meMO 模型 ,但 他 们 提出 的 两 层 建 模 的 思想 和 效应 
量 定义 的 出 发 点 对 调节 效应 拓展 到 更 复杂 的 模型 
是 有 启发 的 。 本 研究 将 借鉴 两 层 模型 定义 的 思路 ， 
将 其 拓展 到 复杂 的 有 中 介 的 调节 模型 的 情景 ， 并 基 
于 此 定义 既 适 用 于 方差 齐 性 条 件 ， 又 适用 于 方差 非 
齐 性 条 件 的 meMO 效应 量 指标 。 

综 上 , 本 研究 主要 解决 以 下 问题 : 首先 从 有 中 
介 的 调节 的 概念 模型 出 发 , 将 Yuan 等 人 (2014) 定 义 
的 两 层 调节 回归 模型 拓展 到 两 层 有 中 介 的 调节 
(2meMO) 模 型 ; 其 次 通过 对 自 变 量 对 与 因 变 量 了 之 
间 关 系 变 异 来 源 的 分 解 ， 给 出 测量 meMO 效应 量 
大 小 的 指标 ; 然后 通过 模拟 研究 评估 模型 参数 以 及 
效应 量 估计 的 精度 ， 并 通过 与 传统 meMO 比较 ， 验 
证 拓展 模型 的 合理 性 ; 同时 采用 实际 数据 案例 和 
Mplus (Muthén & Muthén, 2017) 语 句 介 绍 2meMO 
模型 的 应 用 以 及 效应 量 的 估计 和 解释 ， 以 方便 应 用 
者 使 用 ; 最 后 指出 方法 的 可 拓展 性 以 及 效应 量 指标 
的 优 缺 点 和 实际 应 用 中 应 注意 的 问题 。 


2 有 中 介 的 调节 模型 的 拓展 


虽然 meMO 模型 可 以 分 为 meMO-I 型 和 
meMO-IIT 型 , 但 从 可 解释 性 和 变量 具体 意义 的 角度 ， 
meMO-II 模型 更 有 意义 ， 有 很 多 研究 者 认为 
meMO-I 模型 中 对 应 的 变量 没有 实际 意义 , 应 用 中 
不 应 被 关注 (Edwards & Lambert, 2007; Hayes & 
Preacher, 2013). 因此 ,本文 重点 讨论 更 具 解 释 性 和 
理论 意义 的 meMO-II 模型 (Kwan & Chan, 2018; 
Hayes, 2018)， 其 对 应 的 概念 模型 见 图 1。 图 1 中 如 
RZ T X >Y ARR, H ZEX >Y ERJ 
节 效 应 被 变量 M 所 中 介 , 这 时 就 产生 了 有 中 介 的 
调节 ( 粗 线 条 标 出 的 间接 路 径 )， 变 量 M 称 为 中 介 的 
调节 变量 。 典 型 地 , M 可 能 只 是 部 分 中 介 了 Z 的 调 
TU, 所 以 Z 在 XY 一 了 上 的 直接 路 径 被 保留 。 图 
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的 概念 模型 只 是 各 变量 之 间 关 系 的 理论 表达 ,不 


(后 面 我 们 再 详细 讨论 这 一 问题 )。 2meMO 模型 对 应 


1 
能 直接 进行 分 析 。 只 有 将 其 转换 为 能 够 反应 概念 模 
型 特征 的 统计 模型 ,才能 进行 进一步 的 数据 分 析 。 


ey 


图 1 有 中 介 的 调节 效应 的 概念 模型 ( 粗 体 线 表示 有 中 介 
的 调节 效应 ) 


本 文 借助 于 Yuan 等 人 (2014) 的 思想 , 采用 两 层 
建 模 的 思路 构建 图 1 对 应 的 统计 分 析 模 型 。 图 1 的 
概念 模型 重点 关注 了 对 了 影响 程度 如 何 受 到 调节 变 
量 的 影响 ， 即 对 于 不 同 的 个 体 i, 了 对 了 的 影响 程度 
如 何 随 着 其 在 调节 变量 上 取 值 的 不 同 而 变化 。 第 一 
层 模 型 描述 meMO 分 析 中 自 变量 蕊 对 因 变 量 Y 
响 程 度 的 变化 ; 第 二 层 描 述 Z 和 MM 对 五 和 了 之 间 
关系 的 影响 ， 以 及 调节 变量 Z 和 中 介 的 调节 变量 MM 
之 间 的 关系 。 值 得 注意 的 是 , XJ XA Y ZA 
的 变异 , 在 考虑 了 模型 中 涉及 到 的 调节 变量 Z 和 中 
介 的 调节 变量 M 外 , 可 能 还 会 存在 没有 解释 的 变 
异 ,因此 在 模型 中 增加 蕊 一 了 影响 的 随机 误差 项 。 
2meMO 模型 可 表示 为 : 


Bl: 
¥,=d,+¢X, +6, (2) 

E2: 
d; = yao + YaMj + Ya2Zi t gi (3) 
Ci = Yeo + Y aM; + Yo5Z; + ug (4) 
Mi = uo * uiZi + ui (5) 


方程 (2)~(4) 描 述 了 自 变量 X, 调节 变量 Z 和 中 介 的 
调节 变量 M 与 因 变 量 了 之 间 的 关系 ; 方程 (5) 描 述 
了 调节 变量 Z 和 中 介 的 调节 变量 M 之 间 的 关系 。 
其 中 yzo M ymo 分 别 表示 因 变 量 了 和 中 介 的 调节 变 
fM BIA ARE; Ya Yao Yoo VeVe 和 yw 为 回归 
AM, Bux X, Z 和 oM 是 均值 中 心 化 的 变量 ， 即 
E(Z;) = E(X;) = E(M;) =0 0 RÆ ey, = e; tug; emi 
Mu, 服从 正 态 分 布 ， 且 E(ey,) = E(ew)= E(us;) = 0, 
var(ey;) = Oey var(ey;) = oom ,Var(u,;) = Ow o emi» ey; 和 
ui 是 相互 独立 的 。 值 得 注意 的 是 ,模型 虽然 用 两 层 
的 形式 表达 , 但 是 实际 数据 并 不 是 府 套 的 两 层 数据 ， 
在 参数 估计 上 没 法 区 分 a; 和 wy， 因此 将 其 合并 为 
ex。 男 外 ,在 模型 的 定义 中 ,允许 c; 存 在 随机 变异 ， 
以 此 来 处 理 误 差 方差 随 X; 变化 的 方差 非 齐 性 问题 


的 两 层 统计 模型 可 用 图 2 表示 。 


图 2 拓展 的 有 中 介 的 调节 (2meMO) 分 析 对 应 的 分 层 统 
计 模 型 ( 粗 体 线 表示 有 中 介 的 调节 路 径 ) 


第 一 层 模 型 描述 meMO 重点 关注 和 解释 的 关 
系 是 对 了 的 影响 ， 其 中 截 距 dd, 和 和 斜率 c 为 随机 系 
数 ; 第 二 层 模 型 进一步 解释 Z 和 MM 对 随机 系数 4d; 和 
c; 的 解释 ， 同时 第 二 层 模 型 还 定义 了 Z 和 MM 之 间 的 
KR. 具体 来 讲 , 方程 (3) 描 述 了 模型 中 了 对 回归 
的 截 距 ;受到 的 影响 ( 除 陪 外 直接 指向 因 变 量 Y 
的 变量 ); 方程 (4) 描 述 了 针对 了 的 回归 系数 c 受 到 
的 调节 (指向 了 到 了 回归 线 上 的 变量 )， 其 中 yy, 表示 
XXT YWER, ya 表示 M 对 对 和 了 之 间 关 系 的 
VA, ya 表示 Z 对 XX 和 了 之 间 关 系 的 直接 调节 ; Jr 
程 (5) 描 述 了 M 对 Z 的 回归 。 这 种 采用 两 层 定 义 模 
型 的 方式 突出 了 meMO 模 型 中 解释 的 重点 是 和 了 
之 间 的 关系 ( 层 ID) 如 何 被 调节 (直接 调节 或 间接 调 
节 )， 从 而 可 有 效 区 分 meMO 与 moME 所 关注 问题 
的 差异 。 

下 面 通过 合并 模型 的 表达 进一步 明确 2meMO 
模型 和 传统 meMO 模型 的 联系 和 区 别 。 将 (3) 和 (4) 
代入 到 (2) 中 ,得 到 合并 模型 的 表达 式 为 : 

¥,=d;+¢X;+& = Yao *YaMi * %a2Zi 


tug + (Yoo + YAM; + YoZi tuci) Xi + & 

= Yao * YaMi + yasLi + YeoXi + Ya X iM; 

+ 2X2; t ug X; + ey, (6) 
再 将 等 式 (3) 代 和 人 到 等 式 (6) 中 ， 得 到 : 
Y, = (V40 + Ya uo) + (Yao + YaYuZi; 

* co t YaYuo)X; + (i + Ya Yw0X;Z; 

+ yqremi + Ya X;eyi + Ui Xj + ey, 

7 (Ya, + YaYuo) * Va, + YaYuZi 

Heo + YaYuo) + (Wc t Ya w)Zi 

* (yerem + uci)] Xi + (yarem + ey). (7) 


需要 说 明 的 是 , 相 比 传统 的 meMO BERI, 方程 
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(6) 和 (7) 多 了 wu,X; 这 一 项 ， 此 项 可 以 有 效 区 分 模型 
解释 的 重点 是 有 中 介 的 调节 ,而 非 有 调节 的 中 介 
(我 们 将 在 讨论 部 分 进一步 解释 meMO 和 moME 的 
区 别 )。 男 外 ,合并 后 的 方程 (7) 中 隐 含 着 一 XM 
和 XZ — XM 两 条 路 径 以 及 对 两 条 路 径 系数 的 限定 
( 见 图 2 和 图 3), EIRE X — XM 的 路 径 系数 与 M 
Xpox EIRE >M ) 相 等 , IIA yuo; 限定 
XZ > XM 的 路 径 系数 与 Z 2 M 相等 , WA © 

其 原因 如 下 : 从 方程 (7) 可 看 出 , X 对 了 的 影响 为 
YotYaYuo, BIE X Xt Y WERI yo, DAR X 
通过 XM 对 了 的 间接 影响 yamo, XM XT Y BUS] Ay 
Ya, WARE Xo XM H yuyo, 81x 
XM >Y X vayyo ; MÆ, XZ Xp 了 的 影响 
YatYaYu 包括 XZ XT Y WEIRIN ya, 以 及 XZ 


通过 XM 对 了 的 间接 影响 yywyi ,XM 对 了 的 影响 为 
Ya ， 可 通过 限定 了 XM 为 ym 使 得 
XZ > XM >Y JI yo t Yayu o. 根据 方程 (6) 和 (7)， 
用 于 2meMO 分 析 所 对 应 的 合并 统计 模型 如 图 3 所 
示 。 对 比 图 2 和 图 3, 图 2 的 解释 更 加 直接 。 


ey; 


图 3 拓展 的 有 中 介 的 调节 (2meMO) 分 析 对 应 的 合并 统 
计 模 型 ( 粗 体 线 表 示 有 中 介 的 调节 路 径 ) 


图 3 中 到 了 路 径 系 数 上 的 实心 圆 点 表 对 应 的 
回归 系数 含有 随机 效应 we ， 用 以 描述 了 的 残 差 方 
差 可 以 随 X 的 变化 而 变化 (方差 非 齐 性 )， 传 统 
memo 模型 中 不 包含 这 项 。2meMO 模型 允许 包含 


统 的 meMO 对 应 的 统计 模型 (Grant & Berry, 2011; 
Liu et al, 2012)fH IE, 传统 模型 中 了 一 XM 和 
X 2 XM 的 路 径 自由 估计 ， 而 在 我 们 推导 的 统计 
模型 中 这 两 条 路 径 分 别 限定 为 与 Z 对 M 回归 的 斜 
率 和 截 距 相 等 ， 对 于 模型 的 自由 度 增 加 了 2。 除 了 
随机 效应 wu 外, 虽然 2meMO 对 应 合并 统计 模型 中 
的 其 他 部 分 与 Kwan Fil Chan (2018) 从 另 一 视角 推导 
的 meMO 的 统计 模型 相同 , 但 两 者 建构 和 解释 
meMO 影响 效应 的 思路 存在 本 质 差异 。2meMO 分 
层 统计 模型 (图 2) 与 概念 模型 的 对 应 关系 更 清晰 ， 
模型 中 关于 和 了 之 间 关 系 变异 的 解释 也 更 易于 理 
fife, 可 以 围绕 和 和 了 之 间 关 系 变异 (cj) 的 解释 ,定义 
调节 变量 Z 的 总 调节 效应 、 直接 调节 效应 以 及 间接 
调节 效应 ， 而 且 这 一 定义 效应 量 的 思路 同样 适用 于 
u,; =0 的 模型 。 


3 meMO 效应 量 的 定义 


在 方程 (2)~(5) 所 定义 的 2meMO 模型 中 , 将 (5) 

代入 (4), 可 以 将 对 了 总 的 影响 效应 c 表示 为 : 
€; = Yeo * YaM; + yesLi tug 

= Yeo + Ya wo * YuiZi Femi)? YZi tue 

= (Veo + Ya uo) + Uca * Ya Zi + Yaeui +U (8) 
第 一 项 yo+yaymo 47h X XI Y Koei, FLA y de 
IR X XT Y WARI, yayu KR X it XM 对 了 
的 间接 影响 ,这 一 项 与 调节 变量 Z 的 取 值 无 关 。 第 
二 项 (yi  YaY ui) 表示 调节 变量 Z 对 了 下 了 的 总 调 
节 效 应 , 其 中 y, 是 Z 对 让 了 的 直接 调节 效应 ， 
Vayu Æ Zii M XT x 9 v frr ggg aw Or 
应 概念 模型 (图 2) 中 粗 体 线 表 示 的 Z 通 过 MM 对 XY 和 
了 之 间 回 归 系 数 的 影响 ， 对 应 于 合并 统计 模型 (图 3) 
中 的 XZ > XM Y), yayu BIA meMO 模型 关注 
的 有 中 介 的 调节 , WN memo 。 推 断 统计 中 将 中 介 
调节 效应 的 检验 转换 为 检验 Ho :yowyo = 0 AE TS Lb 
著 的 问题 。 实 际 中 , 经常 通 过 对 调节 变量 Z 取 给 定 
的 值 z (如 高 于 或 低 于 均值 1 个 标准 差 ), 计算 得 到 
Yaym1z， 来 解释 调节 变量 Z 取 不 同 值 时 , Z 通过 中 
介 变 量 M 对 XY 和 了 之 间 关系 的 影响 程度 (简单 meMO)。 
类 似 于 调节 效应 中 简单 效应 差异 的 检验 ， 对 于 给 定 
的 不 同调 节 效 应 值 Z =z MEZ, 9 z,, 可 得 到 中 介 调 


随机 项 wi 不仅 不 再 需要 误差 齐 性 的 假设 , 而 且 更 
符合 实际 情况 。 因 为 除了 模型 中 考虑 到 的 调 方 变量 ， 
和 和 了 之 间 的 关系 还 有 可 能 依赖 于 模型 中 没 包 含 的 
其 他 调节 变量 ,对 了 影响 的 随机 误差 项 也 合理 
的 刻画 了 这 一 点 。 除 了 这 一 随机 路 径 ,该 模型 与 传 


节 效 应 差异 dif emo = Yam (4 — 22), Scb rs ay Xt 
Hy, :dif eyo =0 进行 检验 ,以 回答 调节 变量 Z BUR 
EME z, 和 z, 时 ,中 介 调 节 作用 是 否 存在 差异 。 
对 对 了 影响 程度 的 变异 ， 即 c; 的 方差 为 : 
var(c;) = vat( yeo + Y4M; + Ye2Zi tuu) 
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> (Yaa +Va Ym) var(Zi)+ yoo *o, (9) 
下 面 我 们 基于 c 变异 的 不 同 来 源 来 定义 可 以 
测量 meMO 大 小 的 效应 量 。 第 一 个 效应 量 
(co tarm) var(2) = (0) 
(Yea  YaY a). Vat(Z;) + YaO2y + ne 
为 调节 变量 Z 对 回归 系数 c 总 变异 所 解释 的 比例 ， 
可 以 用 来 测量 站 对 了 影响 的 变异 中 ,调节 变量 Z 所 
能 解释 的 总 变异 所 占 的 比例 。 效 应 量 
yeqvar(Z;) 
(Ye? AVI var(Z;) + 402 + Ow 
则 描述 了 变量 Z 的 直接 调节 效应 对 回归 系数 c 总 变 
异 所 解释 的 比例 , 可 以 用 来 测量 调节 变量 Z 对 
蕊 一 了 路 径直 接 调节 效应 的 大 小 。 而 效应 量 


Dmemo i uo ind 


Puo toi E 


(11) 


Puo dir T 


(ayu) var(2,) (12) 

(Ye + Yar) var(Z;) +AT + Owe 
描述 了 变量 Z 通过 中 介 的 调节 变量 M 的 间接 调节 
效应 (也 叫做 有 中 介 的 调节 ) 对 回归 系数 c 总 变异 所 
解释 的 比例 ,可 以 用 来 衡量 调节 变量 之 通过 中 介 的 
调节 变量 M 对 对 一 了 路 径 系数 的 调节 效应 的 大 小 。 
值得 注意 的 是 ， 基 于 变异 分 解 所 定义 的 效应 量 ， 由 
于 直接 调节 效应 对 应 的 路 径 系 数 y, 与 间接 调节 效 
应 对 应 的 路 径 系 数 y yu 均 与 变量 Z 有 关 , 两 者 之 间 
不 具有 独立 性 。 因 此 , 通过 路 径 系 数 定义 的 变异 解释 
率 不 具有 可 加 性 ， 其 直接 调节 效应 量 与 间接 调节 效 
应 量 的 和 并 不 等 于 总 的 调节 效应 量 ， 即 : piyo tor # 

Puo dir + ÓuO. ind © 

当 方 差 齐 性 假设 被 满足 时 ， 回 归 系 数 c 的 残 差 
随机 部 分 为 零 , 即 w =0， 等 式 (2)~(5) 中 的 模型 被 
还 原 为 传统 的 有 中 介 的 调节 模型 (简称 为 meMO), 

上 面 定义 的 公式 (10)~(12) 被 简化 为 以 下 公式 。 


总 调节 效应 : 
G) 三 (Yez avua var(Z;) (13) 
MO tot 2 2. 2 
(Yer * YaYu). Var(Z;) + yA im 
直接 的 调节 效应 : 
(f) _ y2, vat(Z;) (14) 
MO dir — 
= (ez Faya) Var(Z;) + yom 
中 介 的 调节 效应 : 
N L0) 
Premo MO ind 
(Yau) var(Z;) (15) 


E (Yer + %a¥ur) Var(Z;) + Y402u 
KE, 我 们 使 用 上 标 () 来 区 分 定义 中 不 包含 的 了 对 
7 未 解释 的 随机 残 差 。 


公式 (13)~(15) 同 样 适用 于 2meMO 模型 ， 其 分 
母 表 示 的 意义 是 在 不 考虑 工人 一 了 影响 系数 的 随机 
误差 (1 ) 时 , 调节 变量 Z 对 XY 了 影响 系数 的 调 
节 。 其 包含 可 解释 的 调节 (y。, +a) valz) 和 通 
过 中 介 的 调节 变量 M 间接 未 解释 的 调节 y2o2 两 
部 分 。 在 可 解释 的 调节 中 包含 调节 变量 Z 直接 的 调 
节 和 通过 中 介 的 调节 变量 M 间接 的 调节 。 因 此, 不 论 
是 误差 方差 满足 齐 性 的 条 件 还 是 误差 方差 依赖 于 站 
的 情况 , 我 们 都 给 出 了 测量 有 中 介 调 节 效 应 的 效应 
量 指标 。 公式 (13)~(15) 同 时 适用 于 meMO 和 2meMO 
模型 ， 而 公式 (10)~(12) 仅 适用 于 2meMO 模型 。 

由 于 本 研究 拓展 的 2meMO 模型 随机 部 分 参数 
更 复杂 , 极 大 似 然 估计 方法 会 遇 到 基本 假设 (误差 
方差 齐 性 ) 被 违背 的 问题 。 考 虑 到 贝 叶 斯 佑 计 在 处 
理 复杂 模型 ,尤其 是 售 有 更 复杂 的 随机 变异 模型 时 
更 加 灵活 的 优势 (Wang Preacher, 2015; Muthén & 
Asparouhov，2012; Asparouhov & Muthén, 2020), 
我 们 使 用 贝 叶 斯 方法 估计 2meMO 中 的 参数 及 效应 
量 。 为 了 便于 比较 ， 贝 叶 斯 估计 方法 也 用 于 meMO 
模型 的 参数 估计 。 在 指定 模型 和 先 验 分 布 后 ,采用 
马尔 可 夫 链 蒙特 卡 洛 (Markov Chain Monte Carlo, 
MCMC) Gibbs 采样 获得 模型 参数 的 后 验 分 布 (Gilks 
et al, 1996), XF Gibbs 采样 收敛 性 的 检查 , 使 用 
Gelman-Rubin 势能 缩减 (potential scale reduction, 
PSR) 统 计量 。 根据 Brooks 和 Gelman (1998) 的 研究 ， 
WMR PSR 统计 量 小 于 1.05， 则 认为 参数 估计 达到 收 
敛 。 对 于 收敛 样本 , 我 们 将 参数 的 贝 叶 斯 估计 定义 
为 其 后 验 分 布 的 均值 ( 称 为 后 验 均值 ， 实际 中 也 可 
用 中 位 数 ),， 并 且 使 用 收敛 样本 的 标准 差 获 得 后 验 
标准 差 。 基 于 每 个 参数 后 验 分 布 2.5% 和 97.5% 的 分 
位 数 ， 得 到 95% 可 信 区 间 (Credibility Interval, 
CD(Song & Lee, 2012)。 


4 模拟 研究 


本 研究 通过 模拟 研究 的 方法 ， 比 较 不 同 条 件 下 ， 
meMO 和 2meMO 在 参数 估计 精度 和 效应 量 估计 方 
面 的 表现 。 通 过 meMO 和 2meMO 的 比较 , 验证 方 
差 齐 性 条 件 不 满足 时 ,拓展 的 2meMO 模型 的 必 
要 性 。 

44 模拟 设计 

对 于 方程 (2)~(5), 参考 Wen 等 人 (2010) 以 及 
Liu 等 人 (2020) 的 模拟 条 件 ,参数 真 值 的 设 定 如 下 : 
(1) 调 节 回归 系数 yy 的 大 小 取 3 个 水 平 : 0、0.2 和 
0.4; (2) 随 机 误差 的 方差 ocx 取 3 个 水 平 : 0、0.25 
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和 0.5; (3) 样 本 量 取 4 个 水 平 :100.200.500 和 1000。 
注意 条 件 o2, = 0 时 , 表示 了 的 误差 方差 满足 齐 性 假 
设 , 这 时 对 应 的 真 模型 为 meMO 模型 。 其 他 条 件 的 
设 定 如 下 : KRE yuo 和 yo 设置 为 0; 将 
X 2Y,22Y,M >Y MIZ > M 的 主 效 应 均 设置 为 
0.4, Bl Ya — 427. = Ym = 0.4; Z; Xt X AY ff 
径 的 直接 调节 效应 设置 为 0.2, 即 7 =0.2 。 自 变量 
了 ,调节 变量 Z ,误差 项 ey 和 ey 服从 标准 正 态 分 
布 N0,1)。 共 有 36 FA, 每 个 条 件 重复 500 次 。 
基于 2meMO 模 型 , 使 用 R 软件 (R Core Team, 2016) 
模拟 生成 数据 。 根 据 公 式 (10)~(15) 可 以 得 到 对 应 于 
ya 和 oz 组 合 下 ,效应 量 的 真 值 ,， 结果 见 表 1。 
4.2 ”数据 分 析 和 评价 指标 

对 于 2meMO 模型 和 meMO 模型 ， 使 用 TAGS 
(Plummer, 2015) 和 R 软件 (R Core Team, 2016) 中 的 
程序 包 “R2jags"* 进 行 贝 叶 斯 估计 ， 共 和 迭代 15000 次 ， 
采用 最 后 5000 次 迭代 获得 参数 的 后 验 分 布 ， 即 内 
燃 (burn in) 的 值 设 定 为 10000。 对 于 模型 中 参数 先 
验 分 布 的 选择 ,基于 文献 建议 (e.g., Song & Lee, 
2012; Wang & Preacher，2015)， 对 于 回归 系数 
Y — (uos usYaos'atsYa2» eos es c2) 采用 独立 正 态 
分 布 的 先 验 ， 对 于 误差 方差 - (02,,02,02.) X 
FA Shar B5 33 D AB AY Fe ee, 这 些 先 验 是 回归 分 析 
和 结构 方程 模型 常用 的 半 共 斩 先 验 (Gelman et al., 
2004)。 具 体 而 言 ， 对 于 回归 系数 ， 设 定 先 验 分 布 
为 : 7%, ~N(0,10°); 对 于 误差 方差 , 设 定 先 验 分 布 
Jj: of -IGamma(0.001,0.001) 。 链 条 数 设 为 3， 为 
了 降低 后 验 分 布 间 的 自 相 关 和 节约 存储 空间 ,将 抽 
取 后 验 分 布 的 间隔 设 定 为 10 (Bl thin=10)。 感 兴 
的 参数 包括 中 介 调 节 效 应 yowo 和 中 介 调 节 效 应 的 
差异 dif,,yo， 取 2Z=1 来 评估 的 简单 meMO 效应 
(memo ) 的 估计 , 取 Z=1 和 2Z=-1 来 检验 调节 变量 
Z 取 不 同 值 时 两 个 简单 meMO 效应 的 差异 


(difnemo )o 
对 于 参数 7 uo 和 dif, uuo, 使 用 参数 估计 的 偏 
差 (Bias) 、 均 方 误差 (MSE) 95% FI fci EX [8] AE TER 
检验 力 和 第 一 类 错误 5 个 指标 评估 参数 估计 的 精 
度 。 对 于 效应 量 , 采用 前 3 个 指标 评估 其 估计 精度 。 
同时 采用 meMO 和 2meMO 模型 进行 估计 ， 对 于 
2meMO 模型 ,可 以 得 到 CD OO aroo 
Guo _ tor» PMO_ dir 和 eo 的 估计 结果 ,对 于 meMO 模 
型 ， 可 得 到 前 3 个 效应 量 的 估计 值 。 
偏差 (Bias) 和 均 方 误差 (MSE) 的 定义 如 下 : 
Bias(0) - 0 -0 (16) 
Xn, - 6» 
R 


MSE(Ó) = (17) 


KB, =X} Â / RRR R=500 次 重复 参数 估计 
的 均值 ，0, 为 第 + 次 重复 的 参数 估计 值 。 

95% 可 信和 区间 的 覆盖 率 是 指 R=500 次 重复 中 
Ji 叶 斯 估计 得 到 的 后 验 分 布 中 第 2.5% 和 第 97.5% 
可 信和 区 间 内 包含 参数 真 值 的 次 数 所 占 的 比例 。 检验 
力 是 在 参数 对 应 的 真 值 不 为 0 的 条 件 下 , R=500 次 
重复 中 能 够 检验 出 参数 显著 不 为 零 的 次 数 所 占 的 
比例 。 第 一 类 错误 是 在 参数 对 应 的 真 值 为 0 的 条 件 
F, R=500 次 重复 中 错误 的 检验 出 参数 显著 不 为 0 
的 次 数 所 占 的 比例 。 
43 ”结果 

对 于 每 一 组 模拟 条 件 ， 模 型 中 所 有 参数 对 应 的 
PSR 统计 量 均 小 于 1.05, 说 明 对 于 meMO 和 
2meMO 模型 ， 贝 叶 斯 估计 均 不 存在 不 收敛 的 问题 。 
4.3.1 oo 和 好 wo 的 估计 精度 

X 2 给 出 了 有 中 介 调 节 效 应 估计 值 ovo 的 偏 
差 , 均 方 误差 , 95% 可 信 区 间 覆 盖 率 和 拒绝 率 。 对 于 
JEX, ypeno =0 时 拒绝 率 的 值 代表 第 一 类 错误 
AR, meu 0 时 拒绝 率 的 值 代表 检验 力 。 可 以 看 出 ， 
对 于 参数 估计 的 偏差 和 均 方 误差 两 种 模型 得 到 的 


表 1 模拟 设计 条 件 对 应 的 效应 量 真 值 


Ya 0 0.2 0.4 
c2. 0 0.25 0.5 0 0.25 0.5 0 0.25 0.5 
«2 ai 1 1 1 0.6622 0.6622 0.6622 0.4475 0.4475 0.4475 
«a 1 1 1 0.3378 0.3378 0.3378 0.1381 0.1381 0.1381 
Ao 0 0 0 0.0541 0.0541 0.0541 0.0884 0.0884 0.0884 
Üuo ix 1 0.1379 0.7410 0.6622 0.2128 0.1268 0.4475 0.2402 0.1641 
Puo _di 1 0.1379 0.7410 0.3378 0.1086 0.0647 0.1381 0.0741 0.0507 
Prao 0 0 0 0.0541 0.0174 0.0103 0.0884 0.0474 0.0324 
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表 2 ”两 种 模型 得 到 的 yo NRE, HARE, 95% 可 信 区 间 覆 盖 率 和 拒绝 率 


评价 指标 。 a cl Hfi = = 
100 200 500 1000 100 200 500 1000 
Bias 0 0 0 0.002 0.001 0.000 -0.001 0.002 0.001 0.000 -0.001 
0.25 0 -0.002  —0.001 0.000 0.001 -0.001  —0.001 0.000 0.001 
0.5 0 0.003 0.000 0.001 —0.001 0.003 0.000 0.001 —0.001 
0.2 0 0.08 -0.002  —0.002 0.000 0.001 -0.002  —0.002 0.000 0.001 
0.25 0.08 —0.006 0.000 0.000 0.000  —0.006 0.000 0.000 0.000 
0.5 0.08 —0.002 0.002  —0.002 0.001 —0.001 0.000 — —-0.002 0.002 
0.4 0 0.16 0.000 0.001 0.001 0.000 0.000 0.001 0.001 0.000 
0.25 0.16 -0.004  -0.001 0.002 0.000 | -0.004 -0.001 0.002 0.000 
0.5 0.16 —0.002  —0.004 0.000 0.002  -0.002  —-0.003 0.001 0.002 
MSE 0 0 0 0.002 0.001 0.000 0.000 0.002 0.001 0.000 0.000 
0.25 0 0.003 0.001 0.001 0.000 0.003 0.001 0.001 0.000 
0.5 0 0.004 0.002 0.001 0.000 0.004 0.002 0.001 0.000 
0.2 0 0.08 0.002 0.001 0.000 0.000 0.002 0.001 0.000 0.000 
0.25 0.08 0.004 0.002 0.001 0.000 0.003 0.002 0.001 0.000 
0.5 0.08 0.005 0.002 0.001 0.001 0.005 0.002 0.001 0.000 
0.4 0 0.16 0.004 0.002 0.001 0.000 0.004 0.002 0.001 0.000 
0.25 0.16 0.005 0.002 0.001 0.000 0.005 0.002 0.001 0.000 
0.5 0.16 0.006 0.003 0.001 0.001 0.005 0.003 0.001 0.001 
js an 0 0 0 0.952 0.942 0.942 0.942 0.952 0.956 0.944 0.940 
0.25 0 0.908 0.920 0.918 0.872 0.940 0.956 0.942 0.932 
0.5 0 0.912 0.874 0.852 0.870 0.948 0.952 0.936 0.948 
0.2 0 0.08 0.962 0.950 0.966 0.924 0.968 0.954 0.972 0.926 
0.25 0.08 0.920 0.918 0.922 0.908 0.930 0.944 0.956 0.946 
0.5 0.08 0.894 0.882 0.884 0.882 0.928 0.946 0.936 0.944 
0.4 0 0.16 0.960 0.948 0.958 0.950 0.962 0.950 0.958 0.946 
0.25 0.16 0.935 0.937 0.930 0.933 0.960 0.948 0.942 0.962 
0.5 0.16 0.926 0.928 0.896 0.882 0.946 0.954 0.940 0.940 
拒绝 率 0 0 0 0.048 0.058 0.058 0.058 0.048 0.044 0.056 0.060 
0.25 0 0.092 0.080 0.082 0.129 0.060 0.044 0.058 0.068 
0.5 0 0.088 0.126 0.148 0.130 0.052 0.048 0.064 0.052 
0.2 0 0.08 0.438 0.762 0.998 1.000 0.399 0.760 0.996 1.000 
0.25 0.08 0.390 0.650 0.970 1.000 0.314 0.606 0.950 1.000 
0.5 0.08 0.362 0.600 0.882 0.994 0.269 0.488 0.838 0.990 
0.4 0 0.16 0.938 1.000 1.000 1.000 0.926 1.000 1.000 1.000 
0.25 0.16 0.834 0.996 1.000 1.000 0.808 0.952 1.000 1.000 
0.5 0.16 0.794 0.970 1.000 1.000 0.730 0.964 1.000 1.000 


TE: 覆盖 率 低 于 0.9， 第 一 类 错误 率 高 于 8% 的 值 用 粗 体 字 表 示 。 yowo =0 时 拒绝 率 的 值 代表 第 一 类 错误 率 ，y,,wwo 0 时 拒绝 率 的 值 
代表 检验 力 。 


参数 估计 结果 几乎 没有 差异 (差异 在 小 数 点 后 第 四 2meMO 模型 的 结果 明显 更 有 优势 。 对 于 检验 力 的 
位 )。 在 满足 方差 齐 性 的 条 件 下 , meMO 模型 与 结果 ,meM9O 在 小 样本 时 略 高 于 2meMO, 其 主要 原 
2meMO 模型 得 到 的 参数 估计 结果 也 基本 一 致 。 但 。” 因 可 能 是 因为 这 些 条 件 下 , meMO 模型 过 高 的 第 一 
是 当 方 差 齐 性 条 件 不 满足 时 ， 随 着 cx 方差 增加 ， ”类 错误 导致 检验 力 的 虚假 提高 。 随 着 样本 量 和 调节 
对 于 95% 可 信和 区间 覆盖 率 和 第 一 类 错误 的 结果 ， ”回归 系数 xy 的 增加 , 检验 力 提高 。 
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dif emo F Y aevo 所 得 结论 一 致 。 传 统 meMO 和 
2meMO 模型 的 偏差 都 很 小 。 在 方差 齐 性 的 条 件 下 ， 
两 种 模型 中 Bias 、MSE 、 对 应 的 95% 可 信 区 间 的 覆 
盖 率 和 第 一 类 错误 的 差异 很 小 。 然 而 , 在 误差 方差 
非 齐 性 条 件 下 ，2meMO 对 dif wo 进行 估计 更 加 准 


相对 于 meMO 模型 的 优势 变 得 越 来 越 明显 。 

X 4 给 出 了 2meMO 模型 得 到 的 效应 量 办 wo 
的 估计 值 。 从 结果 可 以 看 出 , 在 所 有 条 件 下 ,MSE 
值 都 很 小 ， 办 so 对 应 的 95% 可 信和 区间 的 覆盖 率 也 
接近 95%。 随 着 o2, MEERE, Breo 估计 的 


wa, 同时 2meMO 对 应 的 95% 可 信和 区 间 的 覆盖 率 更 
接近 95%, 2meMO 方法 在 控制 第 一 类 错误 率 方面 
表现 更 好 ， 随 着 02 值 的 增加 ,meMO 的 第 一 类 错 
误 率 有 所 增 大 。 
4.3.2 ”效应 量 的 估计 结果 

对 于 总 体 中 存在 中 介 调 节 效 应 (ye = 0.2 或 0.4 ) 
时 ,可 以 使 用 meMO 模型 和 2meMO 模型 ， 利 用 公 
式 (13)~(15) 得 到 效应 量 AYO. tor BHO ar 和 Glo 的 
估计 值 ; 当 回 归 系 数 c 对 应 的 误差 方差 cx 2 0 时 ， 
可 以 采用 2meMO 模型 ， 利 用 公式 (10)~(12) 得 到 效 
应 量 fuo. buo ar PE nemo 的 估计 值 。 效 应 量 的 结 
果 表 明 , 不 同 测 量 指标 的 结果 表现 出 相同 的 模式 ， 
PREACH LEAH gU n guo 的 结果 。 

表 3 给 出 了 两 种 模型 得 到 的 G(/) 的 估计 值 。 
样本 量 越 大 ，o2 的 值 越 小 ， 


Ao 的 估计 精度 越 
高 。 在 o% =0 的 情况 下 , meMO 模型 和 2meMO 模 
型 间 的 差异 很 小 。 但 是 ， 当 os x 0 时 , 2meMO 的 结 


果 更 加 精确 , 并 且 随 着 cx 值 的 增加 ，2meMO 模型 


精度 越 来 越 高 。 但 是 , 值得 注意 的 是 ， 当 样本 量 较 
INET, Baa 估计 的 偏差 较 大 ,尤其 是 对 于 较 小 的 
有 中 介 的 调节 效应 而 言 。 因 此 , 在 实践 中 , 为 了 获 
得 对 bromo 更 准确 的 估计 ， 可 能 需要 更 大 的 样本 量 ， 
或 者 类 似 于 传统 回归 分 析 中 对 RT 的 估计 , 需要 对 
其 偏 度 进行 校正 。 


5 应 用 案例 


数据 和 模型 
为 了 便于 2meMO 模型 与 meMO 模型 的 比较 ， 
这 里 我 们 采用 Kwan 和 Chan (2018) 研 究 中 所 用 的 数 
据 (下载 地 址 : http://dx.doi.org/10.1037/met0000160. 
supp)。 其 数据 来 源 于 2012 年 OECD (Organization 
for Economic Cooperation and Development) 组 织 的 
PISA 测试 (the Program for International Student 
Assessment，PISA)， 具 体 介 绍 见 Kwan 和 Chan 
(2018) 研 究 中 的 介绍 。 该 数据 样本 量 为 3047, 分 析 
的 变量 包括 : 因 变 量 数学 学 业 成 绩 (MATH)， 自 变量 


5.1 


表 3 两 种 模型 得 到 的 效应 量 go 估计 的 偏差 ， 均 方 误差 , 95% 可 信 区 间 覆 盖 率 


m e 2 -" moME 2moME 
评价 指标 Yel Oe 真 值 
100 200 500 1000 100 200 500 1000 
Bias 0.2 0 0.054 0.014 0.006 0.002 0.002 0.015 0.006 0.003 0.002 
0.25 0.054 0.015 0.009 0.004 0.001 0.016 0.010 0.004 0.001 
0.5 0.054 0.023 0.013 0.003 0.003 0.022 0.012 0.003 0.003 
0.4 0 0.088 0.003 0.001 0.001 0.001 0.004 0.001 0.001 0.001 
0.25 0.088 0.000 0.001 0.002 —0.001 0.000 0.001 0.002 —0.001 
0.5 0.088 0.001 —0.002 0.002 0.001 0.002 —0.002 0.001 0.001 
MSE 0.2 0 0.054 0.003 0.001 0.001 0.000 0.003 0.001 0.001 0.000 
0.25 0.054 0.003 0.002 0.001 0.000 0.003 0.002 0.001 0.000 
0.5 0.054 0.004 0.002 0.001 0.001 0.003 0.002 0.001 0.000 
0.4 0 0.088 0.002 0.001 0.000 0.000 0.002 0.001 0.000 0.000 
0.25 0.088 0.003 0.001 0.001 0.000 0.003 0.001 0.001 0.000 
0.5 0.088 0.003 0.002 0.001 0.000 0.003 0.001 0.001 0.000 
LE E 0.2 0 0.054 0.962 0.940 0.956 0.926 0.958 0.944 0.954 0.922 
0.25 0.054 0.958 0.926 0.922 0.914 0.970 0.950 0.962 0.948 
0.5 0.054 0.958 0.902 0.882 0.874 0.980 0.954 0.952 0.934 
0.4 0 0.088 0.968 0.952 0.948 0.954 0.966 0.952 0.948 0.950 
0.25 0.088 0.964 0.944 0.948 0.952 0.964 0.944 0.948 0.952 
0.5 0.088 0.932 0.908 0.914 0.874 0.950 0.940 0.956 0.932 


ik: 覆盖 率 低 于 0.9 的 值 


粗 体 字 表示 。 
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评价 指标 YA o. 真 值 100 200 500 1000 
Bias 0 0 0 0.040 0.024 0.013 0.007 
0.25 0 0.030 0.016 0.005 0.002 
0.5 0 0.023 0.009 0.003 0.001 
0.2 0 0.0541 —0.002 —0.005 —0.006 —0.004 
0.25 0.0174 0.021 0.016 0.006 0.002 
0.5 0.0104 0.020 0.010 0.002 0.002 
0.4 0 0.0884 —0.009 —0.008 —0.005 —0.004 
0.25 0.0474 0.013 0.011 0.006 0.002 
0.5 0.0324 0.017 0.007 0.003 0.002 
MSE 0 0 0 0.003 0.001 0.000 0.000 
0.25 0 0.002 0.001 0.000 0.000 
0.5 0 0.001 0.000 0.000 0.000 
0.2 0 0.0541 0.002 0.001 0.000 0.000 
0.25 0.0174 0.001 0.001 0.000 0.000 
0.5 0.0104 0.001 0.001 0.000 0.000 
0.4 0 0.0884 0.002 0.001 0.000 0.000 
0.25 0.0474 0.002 0.001 0.000 0.000 
0.5 0.0324 0.002 0.001 0.000 0.000 
[E ES 0.2 0 0.0541 0.932 0.936 0.948 0.922 
0.25 0.0174 0.978 0.954 0.954 0.936 
0.5 0.0104 0.976 0.964 0.940 0.948 
0.4 0 0.0884 0.946 0.946 0.950 0.932 
0.25 0.0474 0.962 0.946 0.936 0.936 
0.5 0.0324 0.960 0.962 0.932 0.948 


Hb: 对 于 2meMO 模型 ， 


ne 
f 


数学 自我 效能 (MEFF)， 调 节 变 量 社会 经 济 地 位 
(ESCS) 以 及 中 介 的 调节 变量 家 庭 教育 资源 (HRES)。 
模型 假设 : ESCS 可 以 调节 MEFF 5 MATH 之 间 的 
关系 , ESCS 较 低 的 学 生 , MEFF 对 MATH 的 影响 较 
we, 而 ESCS 较 高 的 学 生 ，MEFF 对 MATH 的 影响 
较 弱 ; 进一步 ,ESCS 对 MEFF 与 MATH 之 间 关 系 
的 调节 作用 被 HRES 所 中 介 , ESCS 高 的 学 生 拥 有 较 
多 的 HRES， 所 以 MEFF 对 MATH 的 影响 较 弱 。 对 
应 的 概念 模型 见 图 4, 分 析 的 模型 包括 传统 meMO 
模型 和 2meMO 模型 , 分 析 过 程 中 对 MEFF、ESCS 


图 4 ”家庭 教育 资源 对 社会 经 济 地 位 与 数学 自我 效能 与 
数学 成 绩 之 间 中 介 调 节 效 应 模型 ( 粗 体 线 表 示 有 
中 介 的 调节 效应 ) 


于 方差 估计 大 于 0, 置信 区 间 的 覆盖 率 均 不 包含 0。 


fil HRES 做 均值 中 心 化 处 理 。 为 了 帮助 应 用 者 理解 
和 便于 应 用 者 使 用 分 析 程 序 , 我 们 在 网 络 版 附录 中 
给 出 了 默认 先 验 设 置 时 Mplus 8.3 的 语句 。 需 要 再 
次 强调 的 是 这 里 的 数据 是 单 层 次 结构 ,在 Mplus 中 
我 们 借助 于 多 层 模型 采用 贝 叶 斯 估计 方法 实现 
2meMO 的 参数 估计 ,一 个 个 体 就 是 一 个 cluster。 
5.2 分析 结果 

对 于 meMO 模型 和 2meMO 模型 ,在 默认 先 验 
设置 条 件 下 模型 中 所 有 参数 对 应 的 PSR 统计 量 均 
介 于 1.00 到 1.02 之 间 , 进一步 检查 3 条 不 同 链 的 轨 
迹 图 发 现 不 同 链条 收敛 到 同样 的 后 验 分 布 。 表 6 给 
出 了 两 个 模型 的 参数 估计 结果 、 中 介 的 调节 效应 估 
计 值 和 本 研究 提出 的 效应 量 指标 。2meMO 模型 对 
应 的 DIC (Deviance Information Criteria) 为 78000.4, 
meMO 模型 对 应 的 DIC 为 78069.2， 两 个 模型 对 应 
的 DIC 的 差异 ADIC = 68.8, 2meMO 模型 的 DIC 
较 小 ， 根据 Cain 和 Zhang (2019) 所 建议 的 ADIC > 
7 的 模型 比较 标准 ，2meMO 模型 的 拟 合 优 于 


第 3 期 


刘 红 云 等 有 中 介 的 调节 模型 的 拓展 及 其 效应 量 331 


meMO 模型 。 

AR 5 的 结果 可 以 看 出 ， 两 个 模型 得 到 的 参数 
估计 结果 在 大 多 数 参数 上 相差 不 大 。meMO 和 
2meMO 估计 得 到 的 有 中 介 调 节 效 应 的 后 验 均 值 分 
别 为 -2.019 和 -2.165， 对 应 的 95% CI 分 别 为 
(73.439, —0.614)81(—3.757, -0.673), CI 不 包含 0, 说 
明 有 中 介 的 调节 效应 显著 。2meMO 模型 结果 表明 ， 
第 一 层 数 学 自我 效能 对 数学 成 绩 影响 的 总 变异 的 
后 验 均值 var(c)= 484.239，95% 的 CI 为 (229.158， 
737.546); 数学 自我 效能 对 数学 成 绩 影响 的 变异 


有 中 介 的 调节 效应 ) 的 解释 率 的 后 验 均值 为 1.0%。 
分 母 上 去 掉 误 差 u 的 方差 后 ,得 到 效应 量 
BLO s dU. dir FO Dy 的 后 验 均值 的 估计 值 分 别 
为 0.637, 0.275 和 0.106, 95% CI 分 别 为 (0.296， 
1.000)，(0.000,0.667) 和 (0.000,0.206), 说 明 在 不 考 
虚数 学 自我 效能 一 数学 学 业 成 绩 影 响 系数 的 随机 
误差 (4 ) 时 , 在 社会 经 济 地 位 和 家 庭 教育 资源 联合 
解释 的 数学 自我 效能 对 数学 成 绩 影响 的 变异 中 , 社 
会 经 济 地 位 总 解释 率 的 后 验 均 值 为 63.7%, 社会 经 
济 地 位 直接 调节 效应 的 解释 率 的 后 验 均 值 为 27.5%， 


可 以 被 调节 变量 社会 经 济 地 位 和 中 介 的 调节 变量 
家 庭 教育 资源 联合 解释 部 分 Qo + Yan) var(Z;) + 
Y202, 的 后 验 均 值 为 46.643，95% 的 CI 为 (13.412， 
83.211) 。 由 贝 叶 斯 估计 得 到 的 效应 量 uo ww， 
Guro_ air 和 ero 的 后 验 均值 分 别 为 0.059,0.025 和 
0.010, 95% CI 分 别 为 (0.007, 0.122), (0.000, 0.068) 
和 (0.000，0.023), 说 明 在 数学 成 绩 对 数学 自我 效能 
影响 的 总 变异 中 , 社会 经 济 地 位 调节 效应 的 总 解释 
率 的 后 验 均 值 为 5.996, 社会 经 济 地 位 直接 调节 效 
应 的 解释 率 的 后 验 均 值 为 2.5%， 间 接 调 节 效 应 ( 即 


间接 调节 效应 的 解释 率 的 后 验 均值 为 10.6%。 本 例 
H, fumo SO ig 在 数值 上 差异 较 大 ， 其 主要 原 
因 是 因为 所 用 分 母 的 差异 ， 也 从 另 一 角度 说 明了 
2meMO 模型 中 考虑 大 -> 了 随机 误差 的 意义 。 

本 例 中 , 通过 模型 拟 合 的 比较 我 们 倾向 于 采用 
2meMO, 研究 者 可 以 基于 感 兴趣 的 研究 问题 ， 在 
公式 (10)~(15) 所 定义 的 6 种 不 同 的 效应 量 中 选择 合 
适 的 效应 量 ,解释 调节 变量 对 -> 了 之 间 关 系 变 
异 的 解释 。 但 如 果实 际 中 , 我 们 倾向 于 选择 memo 
模型 ， 则 根据 感 兴趣 的 研究 问题 ,在 公式 (13)~(15) 


表 5 meMO 和 2meMO 模型 参数 估计 结果 


统计 是 meMO 2meMO 
后 验 均值 后 验 标准 差 后 验 均值 后 验 标准 差 

路 径 系数 

ESCS—HRES ( f.) 0.484*** 0.011 0.484*** 0.011 

MEFF—MATH ( 7.) 46.438%** 1.385 48.436*** 1.571 

ESCS—MATH (7,,) 14.621 *** 1.687 14.368*** 1.676 

HRES—MATH ( 74) 3.492* 1.672 3.418* 1.644 

ESCS x MEFF- MATH ( f) -2.911* 1.501 一 3.284* 1.640 

HRES x MEFFMATH ( 7) -4.171** 1.478 —4.471** 1.631 

ESCS x MEFF 一 HRES x MEFF ( £,) 0.484*** 0.011 0.484*** 0.011 
残 差 方差 

o2, 6197.923 160.921 5668.715 189.195 

oy 0.749 0.019 0.749 0.019 

on. 一 一 484.239 130.388 
有 中 介 调 节 meMO 

Pura -2.019 0.718 -2.165 0.791 
效应 量 

Quo ot - = 0.059 0.032 

uo ar = - 0.025 0.022 

bo 一 一 0.010 0.007 

e soe 0.623 0.193 0.637 0.192 

0 a 0.261 0.203 0.275 0.206 

UM 0.110 0.057 0.106 0.056 
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中 选择 合适 的 效应 量 。 另 外 ， 虽 然 从 概念 模型 到 统 
计 模 型 的 转换 过 程 中 ,我 们 强调 模型 中 隐 含 了 的 对 
X + XM fll XZ > XM 两 条 路 径 系 数 的 限定 , 但 在 
实际 中 , 这 两 个 限定 条 件 对 于 参数 估计 结果 的 影响 
很 小 ,， 甚至 可 以 忽略 。 因 此 , 为 了 简化 操作 ， 正 如 
Kwan 和 Chan (2018) 指 出 的 实际 中 也 可 以 忽略 这 两 
个 限定 。 


6 讨论 与 建议 


本 文 基于 单 层 的 数据 拓展 了 meMO 模型 , 采 
用 两 层 模型 的 框架 定义 了 可 处 理 误差 方差 非 齐 性 
的 2meMO 模型 ， 基于 此 模型 将 对 了 影响 系数 变 
异 的 来 源 进 行 分 解 ， 并 由 此 定义 了 可 以 测量 总 调节 
效应 、 直 接 调节 效应 和 有 中 介 的 调节 效应 大 小 的 效 
应 量 指标 。 进 一 步 将 贝 叶 斯 方法 用 于 2meMO 模型 
的 参数 估计 ,并 通过 Monte Carlo 模拟 研究 考查 了 
参数 估计 的 精度 ， 并 与 传统 的 meMO 模型 进行 了 
比较 , 验证 了 2meMO 模型 在 处 理 方差 非 齐 性 条 件 
下 的 优势 。 

根据 模拟 研究 的 结果 , 在 方差 满足 齐 性 的 假设 
时 , meMO fll 2meMO 得 到 的 结果 几乎 没有 差异 , 也 
就 是 说 即便 在 传统 meMO 的 假设 下 ，2meMO 也 能 
得 到 无 偏 的 精度 较 高 的 估计 结果 。 但 是 在 方差 齐 性 
的 条 件 不 满足 时 ， 2meMO 得 到 的 参数 估计 结果 在 
95% 可 信 区 间 覆 盖 率 和 控制 第 一 类 错误 上 的 表现 
优 于 memo 模型 ， 这 一 结果 与 方差 非 齐 性 的 潜 变 
量 交互 作用 模型 (Liu et al.，2020) 和 观测 变量 的 调 
节 效 应 回归 分 析 模 型 得 到 的 结果 一 致 (Yang & Yuan, 
2016). 2meMO 模型 的 拓展 为 调节 效应 模型 中 常见 
的 方差 齐 性 假设 条 件 的 违背 (Aguinis & Pierce, 
1998; Alexander & DeShon, 1994) 提 供 了 解决 思路 。 
基于 贝 叶 斯 估计 方法 的 可 实现 性 以 及 模拟 研究 的 
结果 , 我 们 推荐 研究 者 使 用 2meMO 模型 处 理 方差 
非 齐 性 的 情况 ,在 研究 中 报告 与 感 兴趣 的 研究 问题 
相对 应 的 效应 量 指 标 。 

在 有 中 介 的 调节 效应 的 分 析 中 ， 本文 提出 的 
2meMO 模型 和 对 应 的 效应 量 指标 有 以 下 几 方 面 的 
优点 : (1) 采 用 两 层 模 型 定义 meMO 统计 模型 的 方 
式 , 在 概念 模型 的 理解 和 效应 量 的 定义 上 为 如 何 区 
分 变量 的 角色 提供 了 更 为 一 般 化 的 框架 。 可 以 将 模 
型 关注 的 自 变量 和 因 变 量 关系 的 变化 ,定义 为 层 1 
的 模型 ， 然 后 根据 调节 变量 的 定义 , 在 层 2 模型 中 
解释 调节 变量 对 因 变 量 和 自 变量 之 间 关 系 的 影响 ， 
以 及 调节 变量 和 中 介 的 调节 变量 之 间 的 关系 。 这样 


ar 


分 层次 的 模型 表达 , 合理 且 有 效 地 区 分 了 自 变量 、 
因 变 量 、 调 市 变量 和 中 介 的 调 市 变量 的 不 同 角 色 。 
(2) 文 中 提出 的 效应 量 指标 是 基于 自 变 量 和 因 变 量 
之 间 总 变异 ( 即 var(Q)) 的 分 解 所 得 到 的 ， 这 就 为 紧 
密 围绕 “X 与 了 之 间 关 系 的 调 广 ”这 一 关键 问题 , 定 
义 合适 的 效应 量 提供 了 基础 ,基于 此 所 定义 的 效应 
量 是 围绕 meMO 分 析 所 感 兴趣 的 问题 提出 的 , 更 
具 可 解释 性 ， 同 时 也 更 符合 效应 量 概念 的 核心 内 涵 
(e.g., Kelley & Preacher, 2012)。(3) 基 于 2meMO 提 
出 的 效应 量 指标 ,不 仅 适 用 于 方差 齐 性 而 且 适 用 于 
方差 非 齐 性 的 条 件 ， 有 更 广 的 适用 范围 。 研 究 者 在 
使 用 过 程 中 不 用 担心 由 于 实际 数据 不 满足 方差 齐 
性 而 带 来 的 估计 偏差 检验 力 降低 等 问题 。 
(4)2meMO 模型 的 参数 估计 很 容易 在 Gibbs 取样 的 
框架 下 , 借用 已 有 的 免费 贝 叶 斯 估计 程序 实现 。 如 
2meMO 的 参数 估计 很 容易 通过 常用 的 JAGS 
(Plummer, 2015), WinBugs (Lunn et al.，2012) 和 
Mplus 等 软件 实现 。(5) 更 重要 的 是 ， 本文 所 提出 的 
有 中 介 的 调节 模型 的 效应 量 指标 可 满足 应 用 研究 
者 报告 效应 量 的 需求 ， 弥 补 了 目前 尚 没 有 合适 的 测 
量 meMO 影响 的 效应 量 指标 的 局 限 。 

另外 , 本 研究 定义 模型 的 方法 和 效应 量 的 思路 
具有 折 展 性 。 本 文 提出 的 定义 2meMO 模型 的 思路 ， 
可 以 在 概念 上 为 如 何 有 效 区 分 meMO 和 moME 两 
类 模型 提供 一 些 借鉴 。 从 概念 模型 来 看 , 2meMO 模 
型 和 moME-II 强调 的 重点 不 同 ,两 者 并 不 等 价 ， 
2meMO 模型 是 为 了 解释 调节 效应 是 如 何 通过 中 介 
的 调节 变量 发 挥 作 用 的 ; 而 moME-II 是 为 了 解释 
中 介 效 应 是 如 何 依赖 于 调节 变量 取 值 的 。 然 而 , 传 
统 的 方差 齐 性 假设 下 由 于 两 者 对 应 的 统计 模型 相 
同 ,在 结果 解释 上 很 难 有 效 帮 助 应 用 研究 者 区 分 两 
类 模型 的 本 质 区 别 。 在 两 层 建 模 的 框架 下 ,对 于 
moME 模型 可 以 借鉴 本 文 的 建 模 思 路 , 在 层 1 定义 
与 中 介 效 应 有 关 的 路 径 (a 和 b), 而 在 层 2 解释 中 介 
路 径 中 某 个 阶段 的 效应 被 调节 。 例 如 ,含有 第 二 阶 
段 被 调节 的 中 介 模 型 可 以 定义 为 : 
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M; -a;X;tey 
Y, =d; *b;M; te; X; * e; 


a; — Yao 
d; = Yao + YaiZi + Uai 
b; = Vbo + Y Z; + Upi 


Cj = Yeo 
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上 述 模 型 合并 后 的 表达 式 中 包含 Mu 这 一 项 ， 而 
2meMO 模型 中 包含 的 是 Xu;， 这样 统计 模型 的 定 
义 上 可 将 两 者 加 以 区 分 。 即便 上 述 模 型 中 不 包含 随 
机 误差 项 必 ， 这 一 两 层 模型 建构 的 思路 也 有 利于 
研究 者 区 分 两 类 模型 解释 上 的 区 别 。 未 来 研究 中 可 
以 借鉴 本 文 定义 效应 量 的 思路 ， 基 于 被 调节 的 中 介 
所 关注 的 问题 发 展 合适 的 效应 量 指标 ， 以 解决 有 中 
介 的 调节 模型 与 第 二 阶段 被 调节 的 中 介 模 型 经 常 
出 现 的 结果 解释 上 被 混淆 的 问题 。 

Lachowicz 等 人 (2018) 讨 论 了 合理 的 效应 量 指 
标 应 该 满足 的 特性 (Kelley & Preacher, 2012; 
Preacher & Kelley, 2011; Wen & Fan, 2015)。 本 研究 
所 定义 的 效应 量 指标 具有 以 下 特性 。 QU E AE 
MAG meMO 的 概念 模型 匹配 ,在 变异 分 解 的 框架 
F, 效应 量具 有 可 解释 的 尺度 ,可 以 回答 有 中 介 调 
节 对 总 调节 变异 方差 的 贡献 率 。(2) 效 应 量 指标 不 依 
赖 于 样本 量 ， 且 基于 贝 叶 斯 估计 的 后 验 分 布 可 得 到 
效应 量 的 可 信 区 间 (CD。(3) 风 sor PI duo, BIBAR 
介 于 0 和 1 Zi, KBR, d use o. 
Puo air 和 Gero 的 取 值 也 介 于 0 和 1 之 间 ; 但 在 实 
际 应 用 中 ,可 能 会 出 现 这 四 个 效应 量 的 数值 大 于 1 
(Lutz, 1983; Maassen & Bakker, 2001; Tzelgov & 
Henik, 1991)， 此 时 说 明 Z 对 对- 了 的 直接 调节 (y,) 
和 通过 M 对 和 7 的 间接 调节 (yir ) 方 向 相反 ， 出 
现 了 变量 之 间 的 抑制 效应 (suppression effect) 
虽然 所 提出 来 的 效应 量具 有 很 多 优点 , 但 在 应 
用 时 也 应 关注 其 局 限 性 。 首 先 ， 类 似 于 方差 分 析 中 
的 偏 们 ,var(c) 是 定义 效应 量 的 基础 (分 母 )， 即 效应 
量 的 解释 都 是 以 var(ci) 作 为 比较 基准 的 ， 因此 这 些 
效应 量 本 质 上 是 一 个 相对 的 效应 量 指标 。 因 此 , 在 
实际 应 用 中 , 在 报告 效应 量 的 时 候 应 同时 报告 分 母 
的 大 小 ,以 避免 过 度 夸 大 meMO 效应 ， 同 时 也 有 助 
于 了 解 在 X 对 了 的 影响 中 , var(c) 的 大 小 是 否 有 实际 
的 意义 。 研 究 者 应 该 在 关注 var(c) 本 喘 大 小 的 前 提 
下 ,进一步 解释 meMO 效应 的 大 小 。 另 外 , 我 们 建 
议 在 报告 效应 量 的 时 候 ， 同 时 报告 本 研究 提出 的 多 
个 效应 量 指标 ， 为 全 面 了 解 变 量 的 影响 关系 和 强度 
提供 更 多 的 补充 信息 。 同 时 ， 应 该 注意 基于 变异 分 
解 所 定义 的 直接 调节 效应 和 间接 调解 效应 的 效应 
量 指标 不 具有 可 加 性 。 

大 量 的 研究 证 明 不 同 先 验 分 布 会 影响 参数 估 
计 的 精度 (Yuan & MacKinnon，2009)， 因 此 ,实际 
中 如 何 选 择 合适 的 先 验 分 布 是 研究 者 应 该 注意 的 
问题 。 本 研究 选取 了 最 常用 也 是 研究 者 广泛 推荐 的 


无 信息 先 验 (Browne & Draper, 2006)， 结 果 证 明 即 
使 是 无 信息 的 先 验 也 能 够 得 到 准确 的 参数 估计 结 
R, 但 是 这 一 结论 在 实际 应 用 中 并 不 一 定 具 有 普 适 
性 。 已 有 研究 表明 合适 的 有 信息 先 验 可 以 使 得 参数 
估计 精度 更 高 ， Zondervan-Zwijnenburg 等 人 (2017) 
提供 了 如 何 收 集 先 验 信息 的 指导 。 如 果 研 究 者 可 以 
得 到 合适 的 有 信息 的 先 验 , 我 们 也 鼓励 采用 有 信息 
的 先 验 以 提高 参数 估计 的 精度 。 然 而 , 为 了 避免 由 
于 特定 先 验 设置 带 来 的 主观 性 ,应当 考虑 采用 敏感 
性 分 析 进 一 步 确认 贝 叶 斯 佑 计 结 果 并 不 严重 依赖 
于 某 个 特定 的 先 验 分 布 (van de Schoot et al., 2017; 
Zondervan-Zwijnenburg et al., 2017)。 

本 人 研究 只 考虑 了 含有 一 个 调节 变量 Z 和 一 个 中 
介 的 调节 变量 M 的 情况 ,对 于 方差 非 齐 性 也 只 考 
虑 了 结果 变量 了 的 方差 非 齐 性 ,未 来 研究 中 可 以 考 
虑 更 为 复杂 的 meMO 模型 和 更 复杂 的 方差 非 齐 性 
的 情况 。 男 外 ,对 于 本 研究 提出 的 效应 量 指标 ， 如 
何 结合 效应 量 的 特点 , 给 出 效应 量 “ 小 “中 ”、“ 大 ” 
的 判断 标准 ， 也 是 未 来 值得 进一步 研究 的 问题 。 最 
Ja, 模型 中 被 忽略 的 变量 、 变 量 的 测量 误差 和 其 他 
影响 回归 系数 估计 精度 的 因素 也 会 对 效应 量 的 佑 
计 精 度 产 生 影响 , 未 来 研究 需要 对 这 些 影 响 因 素 进 
行 更 为 系统 和 深入 的 探讨 。 
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Abstract 

Mediation and moderation analyses are commonly used methods for studying the relationship between an 
independent variable (X) and a dependent variable (Y) in conducting empirical research. To better understand the 
relationships among variables, there is an increasing demand for a more general theoretical framework that 
combines moderation and mediation analyses. Recently, statistical analysis of mediated moderation (meMO) 
effects has become a powerful tool for scientists to investigate complex processes. However, the traditional 
meMO model is formulated based on the homoscedasticity assumption, which is most likely to be violated when 
moderation effects exist. In addition, routinely reporting effect sizes has been recommended as the primary 
solution to the issue of overemphasis on significance testing. Appropriate effect sizes (ES) for measuring meMO 
effects are very important in reporting and interpreting inferential results. However, there does not exist an 
effective measure that allows us to answer the question regarding the extent to which a variable Z moderates the 
effect of X on Y via the mediator variable (M) in the meMO model. 

The article is organized as follows. First, the two-level moderated regression model proposed by Yuan, 
Cheng, & Maxwell (2014) was extended to a two-level mediated moderation (2meMO) model with single-level 
data, the statistical path diagram was structured according to the conceptual model and the equations. Second, 
several effect sizes were developed for the 2meMO effect by decomposing the total variance of the moderation 
effect. Third, to estimate the parameters of the 2meMO model and the ES measures of the meMO effects, we 
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developed a Bayesian estimation method to estimate the parameters of the 2meMO model. Fourth, a Monte 
Carlo simulation study was conducted to evaluate the performance of the 2meMO model and the proposed ES 
measures against those with the meMO model. Finally, we illustrate the application of the new model and 
measures with a real data example. 

The simulation results indicate that the size of bias and MSE for parameter estimates are small under both 
meMO and 2meMO models whether the homoscedasticity assumption hold or not. The results of the coverage 
rate of the 95% CI for dif,,,; following 2meMO is comparable to those following meMO when the variance 


of moderation error is zero, which is the assumption the meMO model is based. However, when the 


moderation-error variance is nonzero, 2meMO yields more accurate estimates for dif „emo. than meMO does, 


the advantages of 2meMO over meMO become more obvious as the moderation-error variance increases. The 
results of Type I error rate indicate that 2meMO controls Type I error rather well, and the rates are close to 0.05 
or below 0.05 under all the conditions. However, the Type I error rates of meMO tend to be higher than 0.05 
when the moderation-error variance is nonzero. The power rates following the meMO and 2meMO models are 
comparable for the medium or large sample size, or when there is a large difference in meMO effects. While the 
value of power following 2meMO is slightly lower than that following meMO at small sample se, this result is 
mostly due to the inflated Type I error rate of meMO, and larger sample sizes and the smaller moderation-error 


variances correspond to more accurate estimates of Us The results also indicate that, when the 


homoscedasticity assumption of the meMO model is satisfied, the effect size estimates following the two models 
are about the same. However, when the moderation-error variance is not zero, the results following 2meMO are 
more accurate than those following meMO. 

In summary, the article developed a 2meMO model with single-level data and proposed several measures to 
evaluate the size of the meMO effect explained by moderator variables in total, directly, or indirectly. The 
performance of the 2meMO model is compared against that of the traditional meMO model via Monte Carlo 
simulations. Results indicate that, when the assumption of homoscedasticity holds, 2meMO yields comparable 
results with those under meMO. When the homoscedasticity assumption is violated, estimates under 2meMO are 
more accurate than those under meMO. More importantly, the measures of the size of the meMO effect proposed 
in this article can be used as a supplement to the test of meMO effects and will meet the needs for reporting ES 
in practice. Consequently, the 2meMO model is recommended for the analysis of mediated moderation, and the 
effect sizes (ESs) for the interpretation of the effect according to the questions of interest are better reported. 

Key words mediated moderation, heteroscedasticity, effect size, Bayesian estimation 
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附录 : meMO 和 2meMO 模型 贝 叶 斯 
估计 的 Mplus 语句 


1 meMO 模型 
TITLE: mediated moderation model 
DATA: FILE = memo.csv; 
VARIABLE: NAMES ARE ID MATH ESCS HRES MEFF 
CESCS CHRES CMEFF XCESCS XCHRES; 
! CESCS CHRES CMEFF 为 均值 中 心 化 的 变量 ; 
! XCESCS- CMEFF*CESCS, XCHRES= CMEFF* CHRES; 
USEVARIABLES = MATH CESCS CHRES CMEFF XCESCS 
XCHRES; 
ANALYSIS: ESTIMATOR IS BAYES; POINT =MEAN; 
CHAINS =3; PROCESSORS =3; 
THIN = 10; BITERATIONS=(10000); 
BCONVERGENCE = 0.025; 
MODEL: 
MATH ON CMEFF (rc0); 
MATH ON CHRES (rdy1); 
MATH ON CESCS (rdy2); 
CHRES ON CESCS (rml); 
MATH ON XCHRES (rcl); 
MATH ON XCESCS (rc2); 
XCHRES ON CMEFF (rm0); 
XCHRES ON XCESCS (rm1); 
[CHRES MATH] (rm0 ry0); 
CHRES MATH (sig em, sig ey); 
CMEFF WITH CESCS; 
CMEFF WITH XCESCS; 
CESCS WITH XCESCS; 


MODEL CONSTRAINT: 
NEW(var z vl v2 v3 v5 phi fl phi f2 phi f3 ); 
var z-0.935; !z 的 方差 ; 
vl = (re2+rel*rm1)*(re2+1rel*rm1)*var_z; 
v2 7 rc2*rc2*var z; 
v3 = rc1*rml1*rc1*rm]*var z; 
v5 =vl4rel*rel*sig_em; 
phi fl = vl/v5; 
phi f2 = v2/v5; 
phi f3 = v3/v5; 
MODEL CONSTRAINT: 
NEW(memo memo! memo2 dif); 


memo = rcl *rml; 


memol = rel*rm1*0.967; — !z 的 标准 差 为 0.967; 
memo2 = rcI*rm1*(-0.967); 
dif = memol-memo2; 
PLOT: TYPE = PLOT2; 
OUTPUT: CINTERVAL(hpd) TECHS$; 
2 2meMO 模型 的 Mplus 语句 
TITLE: two-level mediated moderation model with single-level 
data (2moME) 
DATA: FILE = memo.csv; 
VARIABLE: NAMES ARE ID MATH ESCS HRES MEFF 
CESCS CHRES CMEFF XCESCS XCHRES; 
USEVARIABLES = MATH CESCS CHRES CMEFF XCESCS 
XCHRES; 
CLUSTER = ID; ! 对 于 单 层 数据 每 一 个 组 只 有 一 名 被 试 
WITHIN = MATH CESCS CHRES CMEFF XCESCS 
XCHRES;  ! 单 层 数据 中 指定 所 有 变量 为 组 内 变量 


ANALYSIS: TYPE IS TWOLEVEL RANDOM; 


ESTIMATOR IS BAYES; 
POINT =MEAN; 

CHAINS =3; 
PROCESSORS =3; 
THIN = 10; 
BITERATIONS=(10000); 
BCONVERGENCE = 0.025; 

MODEL: 

%WITHIN% 

c | MATH ON CMEFF; 
MATH ON CHRES (rdy1); 
MATH ON CESCS (rdy2); 
CHRES ON CESCS (rml); 
MATH ON XCHRES (rcl); 
MATH ON XCESCS (rc2); 
XCHRES ON CMEFF (rm0); 
XCHRES ON XCESCS (rm1); 

[CHRES] (rm0); 

[MATH] (ry0); 
CHRES (sig em); 
MATH; 

CMEFF WITH CESCS; 
CMEFF WITH XCESCS; 
CESCS WITH XCESCS; 


%BETWEEN% 
[c] (re0); 
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c (sig2 uc); 


MODEL CONSTRAINT: 

NEW(var z vl v2 v3 v4 v5 phi 1 phi 2 phi 3 phi fl phi f2 
phi f3 ); 

var z =0.935; ! the variance of z; 

vl = (rc2*rc1*rm1)*(rc2*rcl *rm1)*var z; 

v2 = rc2*rc2*var z; 

v3 =rel*rm1*rel*rm1*var_z; 


v4 = v1-rcl*rcl*sig em-sig2 uc; 


v5 = vl4rel*rel*sig_em; 
phi 1 = vl/v4; 
phi 2 = v2/v4; 


phi 3 = v3/v4; 


phi fl = v1/v5; 
phi f2 = v2/v5; 
phi f3 = v3/v5; 


MODEL CONSTRAINT: 
NEW(memo memol memo? dif); 
memo = rel*rm1; 
memol = rcl*rm1*0.967; 
memo2 = rel *rm1*(-0.967); 
dif = memol-memo2; 

PLOT: TYPE = PLOT2; 

OUTPUT: 

CINTERVAL(hpd); 


TECHS; 


